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Abstract 

Some years ago, Luck proposed a relevance criterion for the effect of aperiodic 
disorder on the critical behaviour of ferromagnetic Ising systems. In this article, we 
show how Luck's criterion can be derived within an exact renormalisation scheme for 
Ising quantum chains with coupling constants modulated according to substitution 
rules. Luck's conjectures for this case are confirmed and refined. Among other 
outcomes, we give an exact formula for the correlation length critical exponent 
for arbitrary two-letter substitution sequences with marginal fluctuations of the 
coupling constants. 



1 Introduction 

The influence of aperiodic (dis-)order on the thermodynamic properties and the critical 
behaviour of spin systems has been an active research area, particularly since the discovery 
of quasicrystals in 1984. In view of the many articles existing, we cannot give proper 
credit to all contributors; instead, we would like to refer the reader to a recent review 



l~3|j which contains a rather complete bibliography on this subject. For ferromagnetic 
Ising systems, heuristic scaling arguments, put forward by Luck flTSfl , lead to qualitative 
predictions on the relevance of disorder implemented through an aperiodic variation of the 
coupling constants which go beyond the original relevance-irrelevance criterion of Harris 
For the Ising quantum chain (IQC), or the equivalent 2d classical Ising model [O 



with layered disorder, this relevance criterion reads as follows: The critical behaviour 
is of Onsager-type as long as the fluctuations in the sequence of coupling constants are 
bounded, whereas for unbounded fluctuations the critical singularities resemble those of 
the randomly disordered case. This behaviour had first been conjectured by Tracy 
on ths basis of his results for two different examples of three-letter substitution rules. Of 



f 



particular interest is the marginal case with logarithmically (in the system size) diverging 
fluctuations. Here, the critical exponents are predicted to depend continuously on the 
coupling constants 



Meanwhile, some exact results were obtained for substitution chains with marginal 
fluctuations of the couplings. For a number of substitution chains, the surface magnetisa- 
tion critical exponent could be calculated []26], |6], [7|, [|; and, recently, the correlation length 
critical exponent was derived for two special substitution rules using an exact renormali- 



sation approach flq| . All these results are in accordance with Luck's predictions. 

The real space renormalisation method we use in this article is based on the decimation 
process introduced in Ref. fllSl . Our setup is quite different, though; and in particular we 



are going to show how arbitrary substitution rules can be analysed in terms of S'-matrices 
and their star-product formalism. In this way, the dependence of the critical behaviour on 
the nature of the fluctuations can be derived explicitly. The analysis is carried out in detail 
for two-letter substitution chains, with the result that Luck's relevance criterion emerges 
from the renormalisation relations in a natural way. An extended renormalisation scheme 
for the most general case of n-letter substitution rules is presented in an Appendix and 
discussed for certain examples. For infinite classes of substitution chains with marginal 
fluctuations, including the two-letter chains, an analytic expression for the correlation 
length critical exponent is given. 

This article is organised as follows. In Sec. 2, we give a brief outline of the model, 
including some introductory remarks on substitution rules and on properties of the cor- 
responding substitution sequences. Then, in Sec. 3, we introduce the renormalisation 
transformation for a restricted class of n-letter substitution chains, comprising all two- 
letter substitution chains. The dependence of the renormalisation flow on the fluctuations 
of the couplings is shown and the consequences are discussed in detail. Subsequently, in 
Sec. 4, we summarise the implications of our results for the critical properties of the IQC 
for the three different types of fluctuations in the sequence of the coupling constants. Our 
conclusions and an outlook on future developments are presented in Sec. 5. The paper 
contains two Appendices related to the content of Sec. 3: In Appendix A, the solution 
of the eigenvalue problem for the case of a two-letter substitution sequence is derived; 
and finally, generalising the discussion of Sec. 3, the renormalisation transformation for 
arbitrary n-letter substitution rules is explained in Appendix B. 

2 Ising quantum chain and substitution sequences 

The IQC in a transverse magnetic field is defined by the Hamiltonian 



/ N N \ 

2 Es^-i-EVi (2-i) 

\j=l j=l / 
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acting on the tensor-product space (£)f =i C 2 = C 2N . Here, the hj denote the transverse 



x,z 



fields at the sites. The Ej are site-dependent coupling constants, and the operators o •■ 
denote Pauli's matrices acting on the jth site. In what follows, we consider closed chains 
with periodic boundary conditions which are defined by u% +1 = of; various other kinds 
of boundary conditions can be chosen without changing the key results. 

For a general set of coupling constants and fields hj, the IQC can be written as a 
free-fermion model via a Jordan- Wigner transformation and can then be diagonalised by 



a Bogoljubov-Valatin transformation, resulting in [|16 

N 



H N = J2 A Mv q ~l) + C (2.2) 



where 77L 77 are N fermionic creation and annihilation operators and C is some constant. 
The dimensionless excitation energies A q (which can be chosen to be positive and ordered, 
i.e., A^y > . . . > A 2 > A 1 > 0, while properly adjusting the constant C) satisfy the linear 
difference equations 

A q %(j) = - £j _^ q {j-l)-h^ q {3). (2.3) 

Normally, one proceeds by eliminating either ^ or $. In this way, the problem of diago- 
nalising the 2 N x 2^ matrix H N can be reduced to the diagonalisation of an N x N matrix. 
From the set of N eigenvalues (fermion frequencies) A q of the latter, the total spectrum 
of the complete Hamiltonian H N is recovered by considering the sums of the elements of 
any of its 2^ subsets. For periodic boundary conditions, the situation is in fact a little 
more complicated; in order to avoid the appearance of a non-local number operator in the 
fermionic Hamiltonian one has to consider so-called mixed sector Hamiltonians. These 
are obtained by combining the even sector (with respect to the operator Ylj=i a j 
commutes with H N ) of the Hamiltonian ( |2.1|) with periodic boundary conditions with the 
odd sector of the same Hamiltonian with antiperiodic boundary conditions (defined by 
cr^ +1 = —of), and vice versa [|, pJ)| . 

The IQC is equivalent to the 2d classical Ising model with a layered modulation of 
the interactions JT7J. The correspondence shows up through an extremely anisotropic 
limit (also called the r-continuum limit |T0[ ) of the coupling constants of the latter, which 
is believed not to alter the critical behaviour. Explicitly, one establishes the following 
relations: The temperature T of the classical model corresponds to the transverse field 
h of the quantum chain, the free energy translates into the ground-state energy and the 
counterpart of the correlation length £y parallel to the layers is the inverse mass gap 
{Ei — -Eo) 1 of the IQC. In our setting, the mass gap is just given by the modulus of the 
smallest fermionic excitation Ai of ( |2.2j ). According to finite-size scaling, we thus expect 
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a scaling behaviour of Ai at the critical point as 



A 1 (iV) = m^/^k^mN) =: rrfk^mN) , (2.4) 

provided the correlation length scales with a power law at all. Hereby, the scaling expo- 
nent z translates into the correlation length critical exponent v\\ parallel to the layers in 
terminology of the 2d statistical system, while v±_ retains its unperturbed Onsager value 
u ± = l. 

The Hamiltonian (ETTj) is critical (in the sense that the energy gap in the excitation 
spectrum vanishes, i.e., A x = 0) in the thermodynamic limit if the geometric mean (rather 
than the arithmetic mean) of the couplings, normalised by the variables of the transverse 
field, is equal to one, i.e., 

l/N 

lim ^■••^ = 1, (2.5) 

compare Ref. [pOfl . The repetitive structure of the substitution sequences we are going 
to consider guarantees that the mean is well-defined [[H], T§\. Let us now assume that 



the coupling constants Ej (j — 1,2, ... , N) are chosen from a set of n possible values e m 
(m = 1, 2, . . . , n), the value e m occuring with a frequency p m (where Ylm=i p m = 1). For 
simplicity, we take the field to be constant, i.e., hj = h = 1. Then, condition (|2.5| ) can 
be fulfilled by choosing n positive real numbers r x ,r 2 , . . . ,r n with r x • r 2 • . . . • r n = 1, and 
setting the couplings to 

£ m = r m Pm 5 (2-6) 

where p — p x ■ p 2 • . . . ■ p n . Unless stated otherwise, this parametrisation of the critical 
surface will be used in the sequel. 

Substitution sequences 

Let us now consider coupling constants e - drawn from a set of values e a , where the label 
a G A runs over the letters of a (finite) n-letter alphabet A. We choose the couplings 
according to sequences generated by iterated application of a substitution rule g : a — » w a 
on A. For the introduction of the concept and some notation, consider a simple example 
on a two-letter alphabet A = {a, b}: 

(fc) . a -> ab _ fl k\ {k) _ 1 ± V4fc + 1 

We start with an initial word, e.g. W = a, and define W m+1 = ^ fe ^(>V m ). M g (k) denotes 
the corresponding substitution matrix (whose elements count the number of letters a and 
b in the words w a and u^), which has the eigenvalues A± . The entries of the, statistically 
normalised, (right) eigenvector v PF , corresponding to the Perron- Frobenius eigenvalue 
A PF = , give the frequencies of the letters in the limit word W^. In order to ensure 
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that, starting from any seed, the letter frequencies converge to unique, positive values 
in the limit chain, we restrict ourselves to primitive substitution rules g throughout this 
paper. For the corresponding substitution matrices M g , this means that the entries of 
M e are strictly positive for any integer £ > £ (and a suitably chosen finite £q > 0). This 
is obviously true for the examples (|2.7| ). 

Associated to a word W m of length (i.e., number of letters) N = |W m |, one obtains 
finite (closed) quantum chains with N sites by setting e- = e a whenever the jth letter 
of the word W m is an a, and e - = e b otherwise, the final letter of W m determining the 
coupling between the last and the first spin in the chain. While A PF , the largest eigenvalue 
of the substitution matrix, describes the asymptotic scaling of the chain length with the 
number of iterations m, N ~ (App) m , the fluctuation behaviour of the exchange couplings 
is connected to the second-largest eigenvalue A 2 (where it is the modulus of the eigenvalue 
that matters, hence by "second-largest eigenvalue" we always refer to the eigenvalue with 
the second-largest absolute value). In order to quantify the fluctuations, let us define 

where e is the average coupling in the limit chain. Fluctuations stay bounded as long as 
| As | < 1, but for \X 2 \ > 1 they grow with iterated substitutions as a power of the chain 
length, g w ~ N 13 , governed by the so-called wandering exponent [3 that is given by 

(2.9) 



log A 



PF 



In the limiting case | As | = 1, stays constant when we prolong the chain by complete 
iteration steps. Generically, however, finite fluctuations within the substitutes w a of the 
letters a will increase the total fluctuation in the vicinity of these points by a certain 
amount 5: g~ > g w + 5 for some word W with \W\ — |f?(W)| < A and fixed A. By 
repeated use of this argument, we thus expect fluctuations that diverge logarithmically 
with the length of the chain (since the length grows exponentially with the number of 
iterations), see also Ref. 0. On the other hand, the same argument also shows that, 
in this case, the fluctuation behaviour is not completely determined by the substitution 



matrix, but may also depend on the details of the actual substitution rule |12|]; and for 
special choices the fluctuations may still stay bounded even though the second-largest 
eigenvalue is | Aa | = 1. 

Coming back to our examples fl2.7| ), let us concentrate on the three cases k = 1 (which 
yields the famous Fibonacci sequence), k = 2 (the period- doubling sequence), and k = 3 
(the binary non-Pisot sequence). As the second-largest eigenvalues of the substitution 
matrices M g (k) fulfill |A_ | < 1, |A_ | = 1, and |A_ | > 1, this includes an example of 
bounded fluctuations (k = 1), a marginal case with logarithmically divergent fluctuations 
(k = 2) and an example of strong fluctuations (k = 3), hence all three interesting cases 
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are illustrated [ I2[ . Another substitution chain that is referred to several times in the 
text is the celebrated Thue-Morse sequence 



(TM) 



a 
b 



ab 
ba 



A4" p(TM) 



1 1 
1 1 



.(TM) 



A 



(TM) 



(2.10) 



3 Renormalisation transformation 



Rather than reducing the set of linear equations 
does Ref 



15 



]3[) to N dimensions, we proceed (as 
with a direct treatment as a 2iV-dimensional eigenvalue problem of the 



real symmetric matrix 



n 



/o 

hi 























/l 2 



// 


















JV-1 ^N-l 

e N _, 
h N 



-N 











h 



N 



(3.1) 



Since Ti is symmetric and bipartite[J the eigenvalues occur in real pairs of opposite sign. 
But changing the sign of an eigenvalue simply means interchanging the corresponding 
creation and annihilation operators in Eq. (|2.2|) , thus we can restrict ourselves to the 
positive part of the spectrum. In what follows, we concentrate on the case of a constant 
field hj = h, where we may set \h\ = 1 by an appropriate choice of the energy scale. Of 
course, one cannot possibly expect analytical results for an arbitrary sequence of couplings. 
However, for substitution sequences, more can be said using an exact renormalisation 
scheme, and, for an interesting class of chains with marginal fluctuations, the scaling 
exponent z (|2.4|) can be determined exactly and given in a closed form. 

Consider now sequences of couplings chosen according to a substitution rule 



Q 



(3.2) 



on an n-letter alphabet A = {ai, . . . ,a n }. Here, the Wi are n arbitrary words of finite 
length in the alphabet A. Indeed, for the case of aperiodic two-letter substitution chains, 
this form imposes no restriction at all, since these can all be generated by substitution 
rules of the form (|3.2|) . To see this, note that any two-letter substitution rule g can 
be transformed appropriately by inner automorphisms, i.e., by replacing g g : —>■ 



: It can be represented in the block- form Ti, 



A 
A 



after an appropriate reordering of the basis. 
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w^ 1 g(ai)w [where a natural choice of the word w is the common beginning, if any, of the 
words g(ai) and ^(02)], and/or by considering g 2 instead. Obviously, this does not alter 
the limit chain. Note that this is not possible for certain periodic chains, namely those 
which are generated by a substitution rule that replaces all letters by the same word (the 
period), or by a multiple of it. It is no problem, however, to treat this case separately. 



The eigenvalue problem for the matrix TC ( |3.1| ) can be reformulated using transfer 
matrices in the usual way. Our renormalisation transformation (RT) is, however, most 
conveniently described in terms of the star-product formalism, borrowed from scattering 
theory ]32| . We therefore introduce new S-transfer matrices, called S-matrices from now 
on, by 

(&)-«"CW (3 - 3) 

with 

re — 1 

•Sjijfe = Sj\j+i * Sj +1 u + 2 * ... * <Sk~i\k = H< Se\e+i , (3.4) 

£=j 

where the star-product (*) of two 2x2 matrices is defined as 

6 ~ C ) * ( Pl n ) = ( 6 °)+— L - ( ErPl 1 • (3-5) 
r p J \ o o J \0 o J 1 — ppi \ or or 1 p J 



In our case, according to Eq. (|3.1| ), the elementary S'-matrices are 



Smi = 1 4a ) ' (3 ' 6) 

where A denotes an eigenvalue^] of TC. 

In order to establish the renormalisation scheme, we introduce a number of additional 
parameters for these S'-matrices. For each letter a,i e A, we define two asymmetry- 
coefficients Kq. and a field variable h a . that enter the n 2 different elementary S-matrices 
through^ 

Sa\a =\ £ ^ K h A ~ £ ? ka : I • (3-7) 



*3 



The RT now just reverses the substitution steps. In our framework, this is simply done 
by building ^-products of the S-transfer matrices that correspond to the words Wi given 



2 So far, we do not have to impose any restriction here, thus A may be any eigenvalue of H. Only later, 
after expanding the renormalisation equations in powers of A [Eq. ( |3.15 ) and below], we shall restrict 



ourselves to eigenvalues that vanish in the limit of an infinite system size. 

3 Note that this labeling of the fields differs from that used in Eqs. fl2.ip and ([O]) in that the field 
variables are now labeled in the same way as the coupling constants (i.e., as the bonds) rather then by 
the sites of the chain. 
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by the substitution rule ( |3.2|) . Schematically, this works in the following fashion 



rv rv 



a. w] w\ 



rv 

\ w i\ 

w\ a- 



w A w 



3 3 



a, wt w u ... w 



T Q 

a- 



S, 



a-\a ■ 



T Q 
a 3 



T Q 



a 



k 



S, 



a j\ a k 



where w\ denotes the kth letter of the word tOj. In this way, we obtain renormalised 
S'-matrices 



<S ai \ aj — S a .\ w i * S w i\ w 2 * ■ ■ ■ * Sj Wi \ 



e^/c+A —e„ 1 h a . 



(3.8) 



' |aj ' ^ s^ha e~ l K~ A J ' 

compare Eq. ( [3.6|) . As an additional condition on the RT, and in accordance with their 
initial values, we keep the product of the 2n asymmetry-coefficients K^. fixed: 



iixa = n 



K'.K„ 



1 . (3.9) 

i=l i=l 

Now, the RT for the different parameters can be read off directly from the ^-product 
relations. Explicitly, we obtain the following recursion relations 



= h 



_ U l,( k ) 
"t...k I l..,h I La; 



c-2 _ A \( k ) 

t .„fc 



A k M._ 



K = h^ ] (3.10) 



A 



(i) 



A„ 



A 



(fc+i) 



A „ k _ + 



h\A a k) 



e\-A k A™. 

2 



A {1) - A 



A 



(fc + l) _ A (k) 



A fc 



a^+ 



• 2 -A >,A[ k) 



A. 



A.. 



A 



(Kl+i) . 



A 



(Kl+i) . 
a.i+ ' 



for the field variables h„ and for A„ , that are defined as 



K± ■= < A > 



A ai± := ^A 



Using Eq. (|3.9| ), we obtain the renormalised fermion frequencies A as 

' n 



1 

2 li 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



,i=l 



Note the following properties of our renormalisation approach: 
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• The RT is exact, since no new parameters are introduced in a renormalisation step. 
Altogether, there are 3n parameters for n different couplings: Effectively 2n — 1 
asymmetry-parameters k£., the n fields h a ., and the eigenvalue A. Note that the 
coupling constants are not renormalised in our scheme. 

• The rescaling factor of the system size is just given by the leading eigenvalue A PF of 
the substitution matrix corresponding to the substitution rule under consideration. 

• The critical point of the model, defined by the vanishing of the mass gap (A = 0), 
corresponds to a fixed point of the RT for A, as can easily be seen from Eqs. Q3.ll]) 
and (fmj). 

• The scaling behaviour in the vicinity of the critical point is governed by the leading- 
order term of the RT as A — ► 0. 

We now proceed to establish the RT of the spectrum to leading order for A — ► explicitly. 
This is most easily done by noting that, in the product ( p.8|) , the common prefactors 
(f — ppi)^ 1 = 1 + 0(A 2 ) of the ^-products (|3.5|) gives rise to higher-order corrections 
only; and therefore, for the present purpose, the denominator can be neglected completely. 
Consider first the renormalisation of the field variables. We obtain 

K = h ai + 0(A 2 ) . (3.15) 

Instead of working directly with the RT for the fields, it is advantageous to consider new 
variables 

yi :=log(^) . (3.16) 

To linear order in A, the RT of t/i then yields the simple matrix form 

y = M t g y, (3.17) 

where M l Q is just the transpose of the substitution matrix of q. Now, take a look at the 
initial conditions t/ ' for y which are yf^ = — 21og£ ai . As stated above, the entries of 
the (statistically normalised) Perron- Frobenius eigenvector v PF are the asymptotic letter 
frequencies Pi := p a .. We thus obtain 



y° ■ v PF = -2 log 



ITS 



(3.18) 



for an arbitrary set of couplings fulfilling the criticality condition (|2.5|) . Hence, y^ is noth- 
ing but an arbitrary (real) linear combination of the (right) eigenvectors of M l (hence 
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of the left eigenvectors of M e ) excluding the Perron-Frobenius one. Thus, for a generic 
choice of critical couplings, it scales with the second-largest eigenvalue, ||y 



[k]\ 



I Av 



un- 



der the iteration y^- k+l ^ = jVf*i/ fe l. Translated into the terminology of the renormalisation 
group, we recognise y as a scaling field transforming according to 



V = (ApfA 



(3.19) 



where the corresponding renormalisation group eigenvalue is just the wandering exponent 
j3 defined in Eq. (|2.9j ). This demonstrates explicitly that power law, logarithmic, or 
bounded fluctuations lead to relevant, marginal, or irrelevant scaling fields, respectively. 
We are going to describe below how this affects the critical behaviour. 
We now proceed to the RT for A^. Defining 



Kl h 2 k 

k=l 



k c-2 
t .4 



\Wi\ 

Yl II yr + tun 

k=l l=\ w l - 



\wi\ ^2 



n 



r 2 
1 wf 



we find, by induction, as the linear part of the RT of A ai ± 

n 
3=1 

Setting Mt ■= P±., this can also be written in a matrix form 



(3.20) 
(3.21) 



(3.22) 



A± = M ± A ± . 



(3.23) 



Since all the components of the matrices M ± and of the vectors A± are positive, the 
vectors converge to the Perron-Frobenius eigenvectors of M ± under iteration of the RT. 
Let / u ± denote the Perron-Frobenius eigenvalues of M ± . From Eq. (|3.14j) , we conclude 
that 

A 2 = h+ij-A 2 . (3.24) 

For further convenience, we transform M + (under conservation of the spectrum) accord- 
ing to = T~ 1 M + T with a diagonal transformation matrix re- 
entries, we find, dropping the index T, 



h 2 5;-. As new 



13 



K| k h 2^ 

p Za, = n 72^ + • 

k=l 



e , 
=1 w l 



(3.25) 



This form shows explicitly that, to linear order in A, only the reduced couplings e a Jh ai 
enter in the RT. Finally, since the rescaling factor of the chain length is the largest 
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eigenvalue A PF of the substitution matrix, we formally obtain the scaling exponent z 
through 

\j = XjN~* , ,3, 6 ) 
Let us now look at the different cases of fluctuation behaviour in more detail. 



3.1 Bounded Fluctuations 

Because y is perpendicular to the Perron- Frobenius eigenvector of the substitution matrix, 
and since all other eigenvalues are smaller than one in modulus, we find \\y^\\ — > for 
all substitution chains with bounded fluctuations, hence h 2 a ./e 2 a . — > 1. Evaluating the 
transformation matrices M at this fixed point of the fields, we obtain from Eqs. (|3.21| ) 
and fT25| ) 

M+ = M =M e (3.27) 

and hence 

H + = \T = A PF . (3.28) 

Thus we have z = 1 for substitution chains with bounded fluctuations, and the low-energy 
spectrum at criticality scales as 

Aj = x~. (3.29) 

Hereby, the A-transformation induces a finite renormalisation of the fermion velocity v in 
comparison to the uniform chain 

• = £5 M < 3 - 30 > 

A^O ' 

where is the image of A^ = A after a A;-fold application of the RT. 

In general, all renormalisation steps contribute in Eq. (|3.30|) . However, if all eigenval- 
ues of the substitution matrix but the Perron-Frobenius one are zero, the fixed point of the 
RT is already reached after the first renormalisation step and v is obtained immediately. 
This is shown below for periodic chains, but also allows to calculate the fermion velocity of 
the Thue-Morse and related chains. Using the equidistribution property of quasiperiodic 
cut-and-project chains, like our first example (|2.7|) with k — 1 (the Fibonacci sequence), 
the limit in Eq. ( |3.30|) can be carried through explicitly for these cases, too. Since this 
has already been done in Ref . |T8| , we do not repeat the arguments here, and merely cite 
the result 

l \ 21 °g( r ) ,a ot\ 



for two couplings parametrised as in Eq. ( |2.6|) . 
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Periodic chains 



In this short part, we briefly show how the known results for the fermion velocity v for 
periodic quantum chains can be recovered easily within our renormalisation scheme. By 
a substitution rule of the form 

q : di — > w , (3.32) 

any periodic chain with period m = can be generated. In this simple situation, we can 
take all the fields h ai = h and asymmetry-parameters re+ = k = (k„ ) _1 to be equal. By 
virtue of the criticality condition ( p.5| ), the reduced couplings transform as h 2 /ej — > 1 [see 
Eq. ( |3.15[ )], and we obtain the (fermionic) Hamiltonian of a uniform chain (e = h = e w i) 
after a single renormalisation step. Note that the additional asymmetry-parameters result 
in an irrelevant similarity transformation of the matrix ( |3.1| ) only. The renormalisation 
factor for A can be deduced from Eqs. (|3.21|) and ( |3.25| ), explicitly it reads (using Eq. (|2.5|) , 
and normalising to e = h — 1): 

A 2 



2 — 2 

m v 



e\K 2 





\w\ \w\ 

= £ k ' ' ' £ k+l-l(mod\w\) (3.33) 
k=l 1=1 

with e e := e The critical scaling of the fermion frequencies now results in 

277"}' 

A J = 1 ^v, j«iV, (3.34) 



where v is given in Eq. ( 3.33| ), and where the well-known scaling behaviour of the uniform 



chain is recovered for v — 1. This result was first derived in the context of the 2d classical 
Ising model with layered periodicity in Ref. |L]. 

In general, there is no closed expression for the fermion velocity v. For the special case 
of a uniform chain of e a -couplings with periodic e^-defects on every mth site, we derive 

sinh[log(r)J/m 

Note that v(r) = 1 for m = 1 (simple periodic chain), but also v (r) — > 1 for m — > oo since 
the parametrisation (|2.6|) implies that r — > 1 in this limit. 



3.2 Marginal fluctuations 

We now consider the case where the second-largest eigenvalue(s) of the substitution matrix 
(and of its transpose) have |A| = 1. Let E\ denote the corresponding m-dimensional 



12 



(joined) eigenspace of Mi, where 1 < m < n — 1. By Eqs. ( |3.17| ) and ( |3.18| ), y^ 
converges to the projection of on E\ under iteration of the RT. The critical surface 
being (m — l)-dimensional, we can parametrise y = y^ in the renormalisation limit as 

in 

y = J2^s(r i )v i (3.36) 

i=l 

with m vectors v % that span E\ and coefficients > 0, i — 1, 2, . . . , m. This implies 

hi 



exp( yi ) = I]r? (3-37) 



-2 

' a j i=l 



and, if the vector components v * are chosen as integers, the entries of the matrices iVf ± 
are polynomials with positive coefficients of m parameters rf 1 , . . . , r^ 1 . If all rj are equal 
to 1, we obviously obtain 

M + (r, = l) = M~(ri = l) = M e . (3.38) 

Furthermore, if y^ converges to an eigenvector of M l g to the eigenvalue A2 = 1, M + and 
M~ are related by 

M+(r 1; . . . , r m ) = M-^\ . . . , r" 1 ) . (3.39) 

To obtain the scaling exponent z from Eq. ( [3.26D , we need to know the Perron- Frobenius 
eigenvalues ^ of M ± . There is no explicit solution to that problem in the general case. 
In what follows, we present the complete answers for a number of more special, yet infinite, 
classes of substitution rules. 



Cyclic permutations 

Let w be any word from the alphabet, and Pk(w) its kth cyclic permutation. Consider 
now a substitution rule g of the form (|3.2j ) with 

Wi = \p ki (w)] li . (3.40) 

One special case is to take equal words Wi — w for all i = 1, 2, . . . , n. The corresponding 
substitution matrices indeed have A 2 = 1 as their second-largest eigenvalue, with degener- 
acy (n — 1), since the column vectors of M g — I differ by scalar factors only. Thus the fields 
stay at their original values and are not renormalised at all. To derive ^ as a function 
of the parameters r 1; . . . , r n _ 1; we simply note that the special form of g transforms this 
property to M ± , not just for the point r» = 1, Vi, but independently of the parameters. 
We thus obtain 

H± = tr(M ± ) -n+1 > 1 (3.41) 



13 



and the scaling exponent is given by 



log (1 - n + Eti gU) + log (1 ~ n + Er=i Pg-g.) r , ^ 

21og(A PF ) ' (3 - 42) 



where the P^ a . are given in Eqs. (3.21) and Q3.25| ) as functions of the couplings. Since 



/i 1 * 1 are polynomials of the parameters rf l with only positive coefficients, and /U + ({?".«}) = 
£* - ({rr }), we conclude £ > 1. As a function of the parameters, z is either purely convex, 
with unique minimum z(ri = 1) = 1 (the limit of the uniform chain), or it is constant: 
z = 1. It is interesting to take a closer look at the latter case. Obviously, this means 
that P^: a . = [JWjft, independently of the parameters. As can be shown, for substitution 
rules of the form (|3.2| ) with (|3.40|) , this is possible if and only if w is a word (or a power 
of a word) with all letters appearing precisely once and if the last letter of Wi is a«. 
This means, however, that the resulting chain is periodic and has bounded rather than 
logarithmically diverging fluctuations. 

Two-letter substitution rules 

As stated above, each two-letter substitution chain can be generated by a substitution 
rule of the form fl3.2| ). For any marginal substitution rule, we can assume that the second- 
largest eigenvalue of the substitution matrix is A2 = +1. If necessary, this can be achieved 
by going to g 2 . As a consequence, in Eq. ( |3.17| ) is eigenvector to the eigenvalue A2 = 1 
of M , and as above the fields are already at their fixed points. 

In this case, the entries of the transformation matrices M are polynomials of a 
single parameter r that determines the critical couplings. The crucial point now is to 
show that - like in the case discussed above - the matrices M have an eigenvalue A2 = 1 
independently of r and of the detailed form of the substitution rule. This assertion is 
proved in Appendix A. Then, the Perron-Frobenius eigenvalues of A4" ± are given as above 
( |3.4ip , and we obtain the critical scaling exponent 

rM _ log(P+ + P& - 1) + log(P- + P bb ~ 1) ,„ 
Z[ ) ~ 21^) (3 - 43) 

where A PF is the leading eigenvalue of the substitution matrix. As in the case discussed 
above, z(r) is constant (z = 1) only for the special case of a periodic chain (with period 
ab). Otherwise, z is a convex function of r, with a unique minimum z(l) = 1 where the 
(reduced) couplings are equal. On the other hand, it is easy to see that the periodic chain 
with period ab is the only two-letter substitution chain with bounded letter fluctuations 
that can be generated by a substitution rule of the form ( |3.2|) and that has A2 = 1 as the 
second-largest eigenvalue of the substitution matrix. We thus conclude that z = 1 if and 
only if the fluctuations are bounded, and z > 1 otherwise. 



14 



Examples 



The period- doubling chain, given by the substitution rule (2/7) with k — 2, falls into both 
classes described above. We consider double substitution steps, eliminating blocks of the 
form aba and obtain the scaling exponent 

z(r) = lQg(rl b 3 g +p /3) (3.44) 

This result was reported before in Ref. [[05]. If the length of the chain is increased by 
simple substitution steps, the coefficients Xj in Eq. (|3.26| ) converge separately for an even 
and an odd number of steps. Effectively, this can be seen as a renormalisation of the 
couplings according to r = e a /sb — > r _1 in a single renormalisation step. Hence we obtain 
xf d (r) = x^ir- 1 ). 

As an example that is not contained in the above classes, consider the three-letter 
substitution 

a — > abc 

q: b -> ba (3.45) 
c — > ca 

with eigenvalues of the substitution matrix Ai,3 = 1 ± a/2 and A2 = 1. The eigenvector to 
eigenvalue 1 of the transposed substitution matrix M l Q is (0,t, —t), and with Eq. Q3.16Q 
we conclude that h 2 a — > e\ and hl c —> SbE c . Evaluating the RT of A± at this fixed point, 
we obtain the critical scaling exponent 



lo 



g (l + y/l + eje^ + log (l + y/l + e b /e)j 
2 log(l + y/2) 



(3.46) 



Note that, for e b = e c , the fluctuations of the couplings are indeed bounded, and we obtain 
z — 1 as it should be. In our last example, the couplings follow the Thue-Morse sequence 
on the sites of the chain, rather than on the bonds as in all cases discussed so far. The 
coupling constants for a site-problem are given as a function of the two adjacent bonds 
- and thus can take in general n 2 different values (2 2 = 4 in this case). It can, however, 
easily be shown that this can be reformulated as a bond-problem with a substitution chain 
on an n 2 -letter alphabet [ET], In this case, the resulting four-letter substitution rule 



is 

a — > acdb 



b b f . (3.47) 

c — > coac 

d — > dbac 



The asymptotic letter frequencies are p a = pd = 1/6 and p b = p c = 1/3; and so we have 
EaEd^le 1 = 1 (for h — 1) as the criticality condition (|2.5|) . A possible parametrisation 
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is e a = q 2 , Ed = r 2 , e& = s, and e c = (qrs)^ 1 . The eigenvalues of the substitution 
matrix are 4, 1, 1, and 0, we thus expect marginal fluctuations, whereas the Thue-Morse 
bond-problem has bounded fluctuations. Following the above steps, we obtain a scaling 
exponent 

log [(gr)^ + (qr)-^} 

z = i 77^ ~ ■ 3 - 48 

log(2) 

Note that for e a = Ed, £b = £ c we obtain another coding of the period-doubling chain, 
while the special case e a = £&, £ c = £d leads us back to the Thue-Morse bond-problem 
with bounded fluctuations (with q = r -1 , and thus z — 1). 



3.3 Relevant fluctuations 

For IA2I > 1, fluctuations diverge with a power law. As a consequence, the vector 
finally scales with | A2 1 under the RT. Since the matrix elements of M are proportional to 
exp(yj), at least a part of them eventually diverge like exp(c|A 2 |) [note that the components 
of y have different sign because of Eq. (|3.18|) ]. Estimating the product of the Perron- 
Frobenius eigenvalues ^ by the product of the corresponding traces, we obtain 

A 2 ~ jl + fi- ~ A 2 I A2 I (3.49) 

and thus 

A ~ exp (-cN 13 ) (3.50) 
where the wandering exponent [3 is defined in Eq. ( |2.9| ). 



Two-letter substitutions 

Let us discuss this in more detail for the case of two-letter substitutions. We take A2 > 1 
with corresponding eigenvector v of iVf* and parametrise the critical couplings as 

Kb 

-f£ = r v ^ . (3.51) 

£ a,b 

Since v a and Vj, have different sign, either one of these two is bigger than one for a fixed 
value r ^ 1, and without loss of generality we assume r Va > 1. Note that we have the 
following relation 

Kl ^2 

JJ JA = r (A 2 -iK . (3.52) 
1=1 »a 

We now show that P^j > exp [(A 2 — l)i> a logr]. From Eq. ( |3.21|) , this is immediately 
clear for P~ a , but it also holds for P<+. To see this, consider the products of subsequent 
couplings of the word w a . The product of the entire word is given by Eq. (|3.52| ). If the 
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last letter of w a is an a, this enters P^ a and we are done. On the other hand, if the last 
letter is a b, consider the product of all couplings but the last h^/el, which is even larger. 
Iterating this argument, we finally arrive at a letter a which proves our estimate of P^ a . 
Using the standard eigenvalue formula, we obtain the scaling of ^ under the RT, and 
thereby 

A ~ exp (— c| \ogr\N 13 ) . (3.53) 

Here, c is some constant that may depends on the sign of log(r), but apart from that is 
independent of r and N. This scaling form was also predicted in ]18| . 



4 Scaling behaviour and critical properties 

In this section, we give a short summary of the results of our renormalisation treatment 
and discuss their consequences for the critical properties of the IQC. 



4.1 Results: Bounded fluctuations 

Here, the behaviour of the low-energy part of the the spectrum at criticality is essentially 
the same as for the uniform chain. The correlation length critical exponents stays at 
its Onsager value z = u\\ = Uj_ = 1. For periodic chains (or periodic variations of the 
interactions along the layers of the 2d statistical system), the results from renormalisation 
show that the smallest fermion frequencies are only altered by a common factor v (fermion 
velocity), thus also conformal properties persist and the whole model is described by the 
central charge c = 1/2 conformal field theory of a massless free Majorana fermion: 

A, - f (4-1) 



with a fermion velocity v given by ( p.33|) .The scaling dimensions x- are odd integers for 



periodic boundary conditions (mixed sector Hamiltonian) and odd half-integers for free 
boundaries ||. Numerical results [JT2J confirm this conformal behaviour even for non- 
periodic chains with bounded fluctuations. 



4.2 Results: Logarithmically diverging fluctuations 

For the critical scaling, we can make the following Ansatz: 

A, ={£j {r) F r { JlN ) (4.2) 

The index r here denotes the dependence on the position on the critical surface. The 
predicted scaling behaviour |1£| is confirmed analytically for all marginal two-letter sub- 
stitution chains and some infinite classes of n-letter substitution chains. j z F r (j,N) con- 
verges to the scaling coefficients Xj if we increase the length of the chain iV by (multiple) 



17 



substitution steps, in other words F r (J, logiV) is asymptotically periodic for iV ^> with 
period logA PF , where A PF is the scaling factor of the system size. Eq. (|3.43|) implies 
1 < z < oo, and z = 1 is only obtained in the limit of the uniform chain where the 
fluctuations actually vanish. We further note that, in the explicitly solved cases, z is in- 
variant under inversion of all couplings. Thus this is a symmetry of the lower part of the 
critical spectrum, while it is easy to see that it is not a symmetry of the entire finite-size 
spectrum. As a consequence, z quadratically approaches the isotropic value z = 1 for 
small deviations of the couplings from the uniform strength e a — 1, as was also predicted 
by Luck |pj. 

In the terminology of the 2d classical Ising model, we have u\\ — z > 1 for any 
marginal disorder within layers of couplings. Since the correlation length critical exponent 
perpendicular to the layers remains v± — 1, marginal fluctuations introduce a relevant 
anisotropy of the model, as was noticed in Ref. ||. Due to anisotropic hyperscaling, a 
negative specific heat exponent 

a = 2- v\\-v± = l-z (4.3) 

is expected. 

Within the renormalisation approach, the length of the chain N is increased by (mul- 
tiple) substitution steps. Thus, it is a priori not clear what happens to the scaling di- 
mensions when we increase N arbitrarily. Indeed, numerical studies (for free boundary 
conditions) show that, if we increase N by arbitrary amounts, while z converges, F r (j, N) 
stays finite but will not tend to a unique limit. This is in contrast to the bounded fluctu- 
ation case. In more detail, the following properties of F r (j, N) have been observed in all 
examples (of two-letter substitution chains) studied: 

• Given the couplings, F r (j, N) appears to be bounded for all N and j. 

• For N > j > 0, the effect of j on F r is a mere shift: F r (j, N) ~ F(log(iV/j)). 

• As a consequence of the RT, iV ^> j, F r (j, log N) is asymptotically periodic. 

• The amplitude of F r vanishes in the limit of the uniform chain, r — > 1, as F r 
approaches continuously the constant value 7T. 

These points can immediately be translated into properties of the integrated density of 
states (IDOS), the inverse of which is given by 

<j-\p)=p z F r {-\ogp), 0<p<l. (4.4) 

While the first two points imply the convergence of this quantity, the third point above 
implies the lower part of the IDOS to be self-similar with a discrete scaling symmetry 

ma(uj) = (j{m z uj) = a{uj/v) (4.5) 
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where logm is the asymptotic period of F r . Mind that in general not the spectrum (the 
lowest fermion frequencies), but only the IDOS displays this scaling behaviour for N ^> 0. 

This scaling property has consequences for the possible gap structure of the critical 
spectrum, where the results of the renormalisation approach at least lead to necessary 
conditions for gap-labelling theorems. But since within this class more spectral properties 
can be found analytically using trace maps, we postpone the discussion to a forthcoming 
publication. 



4.3 Results: Relevant fluctuations 

In substitution sequences with strong fluctuations, as in our third example, the binary 
non-Pisot sequence [k = 3 in Eq. Q2.7|) 1, their divergence is described ]18| by the wandering 
exponent (3 (|2.9| ) which determines the exponential closing of the spectral gaps as 



Aj ~ (A(r) + B(r)N~ a ^ exp (-F r {j, N)A(r) (y)^ • ( 4 - 6 ) 
Here, A(r) is the second moment of the distribution of coupling constants along the chain 



1 18] which for the binary non-Pisot sequence is given by 



A(r) = Y2V13-7 |log(r)| « 0.459 |log (r)\ (4.7) 

and compatibility with the uniform case requires A(l) = and ctj(l) = 1. 

From the results of the renormalisation approach, F r (j, logiV) should be asymptoti- 
cally periodic for N ^> j. This is confirmed by numerical investigations. Indeed, the latter 
suggest the very same properties as given above in the case of marginal fluctuations. In 
particular, we observe F r (j, N) = F r (log(N/j)) for j 3> 0. 

Again, this means that the integrated spectral density converges in the thermodynamic 
limit and its lower part shows the following scaling property at the critical point: 

a^™-")) = ma (u) , u < 1 . (4.8) 



5 Concluding remarks 

Using a renormalisation approach, we obtained analytic results for the critical behaviour 
of a class of Ising models intermediate between uniformly (or periodically) ordered and 
randomly disordered systems. The systems under consideration are ferromagnetic Ising 
quantum chains whose interaction constants are modulated according to substitution 



sequences. Luck's criterion [18| for the relevance of aperiodic disorder on the critical be- 



haviour is fully confirmed for these systems, bringing along increasing evidence that the 
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fluctuation of the interactions is the basic concept that demarcates the Onsager universal- 
ity class from models with weaker critical singularities. For two-letter substitution rules, 
this analysis was carried out in detail, showing how the renormalisation flow is deter- 
mined by the nature of the fluctuations. An exact renormalisation formalism for the most 
general case of n-letter substitution rules was also presented and applied to examples. 

A number of quantitative results have been displayed, mainly for the correlation length 
critical exponent (resp. the scaling exponent z of the mass gap), confirming and refining 
Luck's predictions for bounded, marginal, and relevant fluctuations of the coupling con- 
stants. In particular, an exact formula has been obtained for the scaling exponent of 
arbitrary two-letter substitution chains, containing previous results as special cases. In 
this light, the case of two- letter substitution rules appears to be rather well understood 
by now. 

What remains to be done in the n-letter case is a thorough discrimination between 
quantum (substitution) chains with bounded and with marginal fluctuations of the cou- 
plings. This is connected with the problem of determining all substitution rules with 
second-largest eigenvalue of the substitution matrix on the unit circle, which nevertheless 
result in a chain with bounded rather than marginal fluctuations. Also, the investigation 
of ordered cases that do not stem from a substitution rule might add some insight as they 
can possess a higher degree of complexity. Unfortunately, renormalisation techniques will 
probably be unsuitable here. 

Of primary interest are, of course, extensions to higher dimensions. A natural first 
step is to study two-dimensional statistical systems with two-dimensional rather than 
one-dimensional variations of the interactions. As can be seen from some special cases ||, 
analytic results are possible in principle, but they are not generic and solvable cases might 
actually be misleading. Further analysis is necessary, and some results can be expected 
from numerical approaches, especially since structures generated by substitution rules 
are well-suited for (numerical) renormalisation. First results in this direction have been 
obtained recently |23) and will be extended soon. 
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Appendix A: Solution of the eigenvalue problem for 

two-letter substitutions 



We consider two-letter chains generated by substitution rules 

a —> aw a f . s 

Q '- b - b Wb (A - 1} 

which have A2 = 1 as the second-largest eigenvalue of the corresponding substitution 
matrix. As a consequence, one finds the constraints 

_ #a(w a ) _ #a{w b ) _ # b (w a ) _ #fe(w 6 ) 
Pa — — j j — — — j j — ; Pb — — i i — — — i i — > 

\w a \ \Wb\ \w a \ \Wb\ 

on the words w a and w b and the letter frequencies p a , b in the limit word, where # a (w(3) 
denotes the number of letters a contained in the word Wp. We parametrise the criticality 
condition (|2.5f) s p a a e p b h = 1 (setting h = 1) for the corresponding Ising chain by 

e a = r- l/ \ e h = r k l\ (A.3) 

with coprime integers k and I, such that lp a = kp b . Then k + 1 divides \w a \ and \w b \. 

The scaling exponent of the critical spectrum is a simple expression of the leading 
eigenvalues of two matrices 

/ p± p±\ 

M ± = ^ 6 | (A.4) 

V r ab r bb / 



where the entries are defined in Eqs. ( |3.21 ) and ( |3.25 ). They are polynomials in 



and r . In what follows, we prove that M have an eigenvalue A2 = 1, independently 
of the value of r and of the detailed form of the substitution rule g. Since we have 
M + (r) = M~{r~ l ), it suffices to concentrate on M + . 

We define two polynomials P a,b corresponding to the words w a and w b by 

P>, r-\ u) := (P+ - 1) + uP+ , P\r, r~\ u) := P+ + u{P b \ - 1) . (A.5) 

Taking the frequencies p^ b of the letters as given, the possible words w a and w b are 
subject to the constraints ( |A.2|) . We define an elementary polynomial P e (r,r~ 1 ,u) to be 
the polynomial P a ' b corresponding to a word w e of minimal length fulfilling this condition. 
Since k and / are coprime, each w e consists of k a's and I 6's and leads to one of 
different elementary polynomials possible. 

Let in the following w be w a or w b . We now introduce a graphical representation of 
the word w and the accompanying polynomials. For given w, construct iteratively a step 
function g w : [0; \w\] — > Z as follows: 

g w (0) = 
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Note that the criticality conditions imply 
/ = 3 is shown in Figure |l|. For w = w a (w = Wb 



0. An example for the case k = 2, 
P+ a — 1 (resp. P b + a ) is given by the "sum 



A 



lwl= 10 



Figure 1: Step function g w (x) (|A,6|) of the word w = abbaabbabb. The endpoints of an elemen- 
tary polynomial that can be separated are marked by 'A' and 'B'. 

over the upward steps" of g w (proceeding in positive x-direction), while P^ b (resp. P h \ — 1) 
is given by the "sum over the downward steps", each step contributing a term r 9w ^ k \ In 
our example (for w = w a ), P^ a — 1 = r 3 -|_ r 2 -|_ r 5 -|_ r 4 an d P^ b = r + r _1 + r 3 + r + r 2 + 1. 

In a first step, we now show that P a ' b can be expressed as sums over elementary 
polynomials P e . 

Lemma: For an arbitrary word w fulfilling the criticality condition for given p a ^, there 
are integers ctj and elementary polynomials P/, such that 



P a ' 6 (r,r-\u) 



Y,r a *P?{r,r-\u) 



(A.7) 



Proof: Consider the substrings Sj of w of length k + / starting with the letter wK For 
each Sj containing more than k a's, there has to be an Si with less than k a's, because of 
the criticality condition for w. Since the number of a's contained in successive substrings 
Sj and Sj + i differs at most by one, there is at least one substring consisting of k a's and 
/ 6's. In the graphical representation, this appears as a restriction of g w to an interval of 
length k + / with the endpoints of the graph taking the same value. The corresponding 
term within P a,b is obviously some elementary polynomial times a power of r. We may 
now eliminate this substring from w (cut out the interval in the graphical picture) and 
start the argument again. 
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We proceed by showing that the elementary polynomials differ only in a (polynomial) 
factor that is independent of u. 

Proposition: Given k and I, there are polynomials Pk{r, r _1 ) and Pi(r,r~ l ) such that for 
every elementary polynomial P e 

P e (r, r~\ u) = Q e (r, r" 1 ) [P k (r, r" 1 ) + uP^r, r" 1 )] (A.8) 

where Q e (r, r _1 ) is a polynomial that does not depend on u. 

Proof: An elementary polynomial P e takes the following general form 

( ci c 2 i \ 

r ~ ifc + rl ~ jk + • • • + J2 rkl ~ jk ( A - 9 ) 
3=1 j=ci j=c k J 

where < c\ < c 2 < . . . < < I are integers. Defining 

Pk--=Y,* i = — 1 p ' := 5>" = t^' (A - 10) 

j=0 3=1 

the proof then follows by direct calculation. We remark that with this choice of Pk and 
Pi, Q e can also be shown to be a polynomial in r and r _1 . 

Now, we are in the position to complete our argument. From Eqs. (|A.5|) , ( A..7 ) and 
( |A.8|) , we conclude that 

Paa 1 Pk _ Pba r » -. -. \ 

P = P" = P l" ' l A " ii J 

-T ab -H -r, f,fe — 1 

Thus it follows that det (M + — I) = 0, and hence A2 = 1 is an eigenvalue of M + for an 
arbitrary value of r. The leading eigenvalues of M ± are given by Eq. (|3.41 ). 



Appendix B: Renormalisation formalism for general 

n-letter substitution rules 

In the above discussion, the coupling constants were chosen according to a substitution 
rule of the special form ( |3.2|) . While any (infinite) substitution chain with two letters can 
be generated that way, this is no longer the case for three or more different letters (resp. 
couplings). In this Appendix, we present a generalised version of the renormalisation 
formalism to deal with the general case. Consider a substitution rule of the form 

q: Oj — ► Wi = w il a i w i 2 • (B.l) 

The main problem in the general case is to maintain the recursive structure of the renor- 
malisation procedure within the elimination process. By fixing the first letter of the 
substitutes in Eq. ( |3.2| ), the appropriate couplings of the renormalised chain were given. 
This is no longer possible here; however, we can proceed as follows. 
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Fix any letter dj within the word Wi = g(a>i), thus defining two words wn and Wi2 
as shown in Eq. ( |B.1|) . Note that the letter a, appears in the word u>, (at least for 
a suitably chosen power g k of the substitution rule) for primitive substitutions. 

Redefine the substitution rule as a function on pairs of letters in the following way 

g-^> g : ai(aj) -> a i w i2 w jl (a j ) (B.2) 

Obviously, ( |B.2| ) leads to the same limit chain as ( |B.1| ). The substitution matrix 
of g has a dimension of at most n 2 (pairs a^aj that do not appear in the original 
chain need not be included in g). Actually, Mg is just the substitution matrix of 
the corresponding site-problem, where the type of site follows the sequence, and the 
bond (resp. the interaction along it) is a function of its two endpoints 0, ^5|. Note, 
however, that the couplings here are still attached to the bonds. 

In the renormalisation formalism, now define S'-matrices, fields and asymmetry- 
parameters corresponding to pairs of letters 

b <H\a d ~ !, i_ . ■ (B.3) 

V c cij "'didj t a . i% a . a .i\ I 

Renormalisation again means inverting the substitution process by integrating out 
all degrees of freedom attached to the words wn and . We obtain renormalised 
S'-matrices through 

S a \a = S n .\„,i * S,„i i,„2 *...*S \w i2 \, i *...*S iwni, . (B.4) 

The RT (and their linear orders in A) for the (at most) 3n 2 parameters are then 
obtained analogous to Section 3. 



For a classification of the critical behaviour based on the fluctuations we need to know 
about the spectrum of Mg in dependence of the spectrum a e of M e . For an appropriate 
power of g, it is indeed possible to show that contains a e and that all additional 
eigenvalues are either 1 or 0. We only sketch the proof here. Without restriction of the 
general case, we can assume that the first and the last letters of all words Wi remain fixed 
under the substitution; that is wj = [g^Wi)} 1 for the first letter, and equivalently for the 
last letter. This is always fulfilled by a finite power g r of the substitution rule. We assume 
the n words Wi to have k different first letters and j different last letters. Taking a look 
at the associated pair substitution matrix Mg, it is easy to show that 

tr(Mg) = tr(M e ) + n + kj - k - j . (B.5) 
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Since we always have Mg-m = (Mg) m and the numbers k and j remain fixed because of the 
above assumption, the traces of M^™ and M e m differ only by the constant n + kj — k — j 
for arbitrary m. We conclude that cjg containes a e and the all additional eigenvalues are 
fixed points under any power, thus 1 or 0. 

We will, however, not give a detailed discussion of the general case here, but rather 
illustrate the method by applying it to a special example. Consider the substitution rule 
of the "circle sequence" . This example of a quasiperiodic chain (in the sense that it has 



a pure point Fourier spectrum) has been studied numerically in Ref. [18 



a — ► cac — > abcaccacabcac 
q : b — > accac — *> cacabcacabcaccacabcac (B.6) 
c — > abcac — > cacaccacaabcaccacabcac 



with £> 2 being of the form (|B.2|). The eigenvalues of M g 2 are r 6 , 1, and r 6 , where 



r = (l + v / 5)/2 is the golden mean - we thus expect marginal fluctuations. The asymptotic 
frequencies of the letters are p a = 2 — r, p& = r — 3/2, and p c = 1/2, respectively. We 
choose e a = rs~ T , = rs, and e c = r~ 1 s as the parametrisation of the critical couplings 
e p a a el b e p c c = 1. Since only the five pairs ab, ac, be, ca, and cc appear in the chain, the 
dimension of the pair transfer matrix Mf. is just five. The spectrum consists of the 
spectrum a e 2 and of two additional eigenvalues 1. Identifying pairs of couplings with the 
first one, we clearly have p a b + p ac = Pa, Pbc = Pb, and p ca + p cc = p c . Thus the vector 
y of the logarithms of the reduced couplings, y xy = 2\og(h xy /e x ), is sill perpendicular 
to the Perron-Frobenius eigenvector of Mq and converges in the renormalisation limit to 
a linear combination of two eigenvectors v 1,2 corresponding to an eigenvalue 1 of 
y — > log r v 1 +log s v 2 . While both parameters, r and s, enter the transformation matrices 
M ± of A±, s can actually be eliminated by a similarity transformation. Such a behaviour 
is always to be expected if a parameter is connected to an eigenvector (like v 2 in this case) 
that corresponds to an additional eigenvalue 1 of iW| 2 , not contained in a e 2 . The Perron- 
Frobenius eigenvalues of can be determined explicitly, and we obtain for the scaling 
exponent 



log (5 + 2(r + r- 1 ) + 5 + 2y/(r + r" 1 + 2) 2 + (r + r" 1 + 2)) 



(B.7) 



6 log(r) 

Finally, we like to show an example of a three-letter substitution rule 

a — > bca 

q: b -> beb (B.8) 
c — » ca 

with bounded fluctuations, which nevertheless has A2 = 1 as the second-largest eigenvalue 
of the substitution matrix. This is the bond substitution that corresponds to the site- 
problem given by g s : a — >• aba ; b — >• ba. Fluctuations are indeed bounded here, since 
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we just obtain the Fibonacci chain by identifying a and b, but also by identifying a and 
c. Calculating the scaling exponent by the above method, we obtain z — 1, as was to be 
expected. 
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